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5. A Describing Function of the CMG Nonlinearity 

Using the Analytical Torque Equation 

A describing function of the CMG frictional nonlinearity was 
derived earlier using the straight-line approximated input-output 
relation between the frictional torque T^p. and the CMG angular dis- 
placement 

However, it is possible to derive a describing function for the 
CMG frictional torque using the analytical relation between T^p and 6g. 

It has been established that the frictional nonlinearity of the CMG 
can be described by the square-law relation. 




(5-1) 


where 


= "gF 


(5-2) 


TgpQ = saturation level of T^p 


Y = positive constant 


Carrying out the integration on both sides of Eq. (5-1) yields 




-1 


^^"“"gf ■ '•'gfo^ 


0G>_O 


(5-3) 


6g + 


-1 


i'^'''gf ''"gfo^ 


9g10 


(5-4) 


where C-j and are constants of integration, and 
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Then, C-, and C 2 ane given by 
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where 8^^ and T^p. denote the initial values of Gq and T^p, respectively, 
For a sinusoidal input, 0^ is represented by 

6g = A cos cot (5-9) 

It is important to note that for the input of Eq. (5-9) 9^^ = -A when 


G g ^ 0 , an d 0 = A when 8 g ^ 0 . 


Solving for Tgp and Tgp from Eqs. (5-3) and (5-4), respectively, we 


have 
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where 
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With the describing function method^ the frictional torque may be 
approximated by the fundamental component of the Fourier series. The dc 
component is zero, since the input-output relation is symmetrical 
about the zero-torque axis. 

Thus , 

Tqp ^ A.J sin tot + cos wt 


"V^l ^1 - (1>) 
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Substitution of Eqs. (5-10) and (5-11) into Eq. (5-19) gives 
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In arriving at the last expression it is noted that 

I 

C-j > A 

and C-j = -C 2 over their respective ranges of Equation (5-23) is 
simplified further to 
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Now substitution of Eqs. (5-10) and (5-11) into Eq. (5-20) yields 
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Evaluating each of the integrals in the last equation, we have 
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Since C 2 < -A, C 2 is always negative, and ^ the first integral of 
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where the fact that C 2 is negative has been used. Also, tan tt/ 2 is 
taken to be +» since oit/Z expands from 0 to it/2. 

Similarly, the third integral of B.j in Eq. (5-25) is written 
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In arriving at the last equation, we have recognized that > A and 
have used that tan tt/2 = since in this case ojt/2 expands from tt/2 to tt. 


Thus , 
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The describing function in complex form is written as 


N(A) = (5-31) 

where A^ and are given by Eqs. (5-24) and (5-30), respectively. A 
digital computer program for the computation of N(A) and -l/N(A) is listed 
in Table 5-1. The constant A is represented by E in this program. The 
parameters of the nonlinearity are: 

Tgpo= O.l ft-lb 
Y = 1.38 X lO® 

Figure 5-1 shows the magnitude (db) versus phase (degrees) plots of 
-l/N(A) for y = I. 38 x lO^, 0,69 x lO^, I. 38 x lO^, 0.69 x lO^, and 
I. 38 X lO^, as the magnitude of A varies. Note that as A becomes 
large, the magnitude of -l/N(A) approaches infinity and the phase 
approaches -270 degrees. As A decreases, the magnitude of -l/N(A) decreases 
and the phase approaches -180 degrees. In the limit as A-^-O, ~l/N(A) 
goes to -VyTgpQ^' 


Asymptotic Behavior of -l/N(A) for Very Small Values of A 
Figure 5-1 shows that as A approaches zero, the magnitude of 
-l/N(A) in db approaches 20 "^og-jQEl/yTgpQ ] and the phase is -180 
degrees. The asymptotic behavior of -l/N(A) for very small values 
of A is derived here analytically. 
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.100 LST Ca^^TI^SO^JS riESCRIBir^G FUMCTIOri - EXRCT CtlG MDNLIMERRITY 

1 • 0 U 0 L- D M F' L E X V ? H 1 E> 

1 .50 0 RERL^8 P < E 0 > j P I ? RRD ? !□ ? GRMMR :« EST RRT ? E >Rfl ? R «. TSF I ? TGPH j TGFP » C 1 PE 

jRl ?B1 

1 .800 PERL^S RZ 

2.000 P 1=3. 14159 

3.000 RRD=130..--PI 

4 . 000 TD= .1 
4.500 S=10. 

5 . 0 0 0 G R M MR = S ♦ 1 . 3 8 E 5 

6.000 ESTRRr=l .E-13 

7.000 MP=£ 

5.000 MB=1S 

9.000 P':l) = l. 

10.000 PC2>=5. 

10.200 I'/RITECG j 1 00> 

10.4 00 WRITECt. j 1 01 > 

11 .000 DO 1 J=l!.MD 

12.000 DO 1 1=1 tNP 

1 3 . 0 0 0 E=ESTRRT^P < I > ^ < 1 0 . o < J- 1 > :> 

14.0 00 HR=2’.-^GRMMR<>E<>TQ 

15.000 R= < - 1 . .''‘RR > +IISQRT < <. RR-*-RR+ 1 . > > < RR^RR > > 

16.000 TGFI=R^rn 

17.000 TGFM=TGPI 

18.000 TGFP=-TGFI 

1 9 . 0 0 0 C 1 =£- 1 . GRMMR^ < TGFP^ TD > 

2 0 .000 C2= -E~ 1 . C GRMMR^ c TGFM+TO > > 

21 .0 0 U R 1 = ( -4 . ♦TQ/P I > + C 1 . < p I ♦GRMMRoE > > ♦IiLDG (. < < Cl + E > ♦ <' r.£-E > '■ < <’ r l ~E '* ♦>- 

21.10 0 < C:£+E > > > 

2 1 . 5 U 0 RZ=DLC]G t. c. f. L- 1 ^ (' C'2”E •' •■■■' c! C C 1 ~£3' ♦ <! C2+E '• 

21.600 Rl = <--4 .♦TO-'PI ;-->-(RZ--'vPI^GRMNfl^E> > 

22 . 0 UO E 1 = - 1 , < GRMMR^E > > <► C £ . + C2.- DSQRT < C2^C£>-E>E > -C 1 --'DSQRT < P 1 ♦P 1 -E-frF '■ > 
22.100 R1=R1 E ^ 

22.200 I:l=Bl.--E 

23 . 0 0 0 GM=IiC;MPLX C B 1 j -R 1 > 

24 .000 bV=~l ..-•■■GH 

25 .000 G 1 =RERL C GV > 

26 .000 G 2= R I MRG < G V > 

27 . 0 0 0 GMRG=CRES: C GV > 


28.000 GDE=2 0 . ♦HLDG 1 0 '•! GMhG ) 

£9 . 000 GPHRSE=RR£i^RTRM2( G£ *G1 > 

3 0 .000 I F c; GPHRSE . GE . 0 . ) GPHRSE=GPHRSE-36 0 , 

33 . 0 0 0 WR I TE < 6 ? 1 02 E ? T GF I > gPHRSE » GDB ? GMRG 

34.000 1 CDMTINUE 

35.0 00 1 00 FORMRTf;-' CDNT I MQIJS DESCRIBIMG FUNCTIDM FDR CMG MDMLIMERRI 


36 . 0 0 0 1 01 FDRMRT < » SX * E :• 1 1 X ? TGF I ■<1 fiX ' •' PHR’'” E '' 
riE--/> 

37.000 102 FDRMRTC 1P5E14 ,5> 

39.000 STOP 


1 OX j DB ^ > 9X MRGM I TU 











Dividing both sides of Eq. (5-24) by A and expanding the loga- 
rithmic term into a power series, we have 
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Taking the limit on both sides of Eq. (5-32) as A-K), and using the 
fact that 
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Substituting Eq. (5-33) into the last equation we have 
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Dividing both sides of Eq. (5-30) by A and taking the limit as A 
approaches zero, we have 
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2 - 1/2 

Expanding [1 - (A/C-j) ] into a pov/er series, and using only the 
first two terms, Eq. (5-35) becomes 
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- tim (5-36) 

A-^ yC^ 

Thus , 

•ttm ]/N(A) = -Km 5 . - J ...... (5.37) 

As shown in Fig. 5-1, the gain-phase plot of -1/N(A) as A approaches 
zero is a point which lies on the -180 deg line with a magnitude of 
20 ^ ^O-jqCI/yTqpq^J* 
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Asymptoti c Be havior of -1/N(A) For Ve ry Large Values of _A 
For very larfje values of A> the value of C-j becomes 

tim C, = A + (5-38) 
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Similarly we can show that 
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Thus , 
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As shown in Figure 5-1, the gain-phase plots of -1/N(A) approach 


-270 ° as A ^ for all values of y and T 


GFO* 



6. Computer Simulation of the Simplified LSI System with the 
Analytical Torque Expressions 

A computer simulation of the 1ST system is presented here to 
corroborate the results of the describing function analysis of the 
last chapter. Since the describing function analysis has been carried 
out with the analytical torque expressions for the CMG frictional 
nonlinearity, the simulation model of the nonlinearity also has the 
same characteristics. This model of the nonlinearity is implemented by 
using the expressions for T^p in Eqs. (2-38) and (2-39) with initial 

conditions for and T^p being redefined each time a sign change in 

e occurs. 

G 

The simplified LST system is represented by the block diagram of 
Figure 1-7. The linear transfer function which the nonlinear element NL 
sees is given by 
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JgJ^s 

^+J^KjS^+KjHK^S+KjHKq 

Two sets of numerical 

values are 

considered as follows: 


System 1 

System 2 


10^ 

10^ 


2.1 

3.7 


216. 

280. 


9700. 

10000 

H 

600 

200 


1371.02 

3000 

^0 

5758.35 

20000 



The frequency- domain plots of 6(s) for both systems are given 
in Figure 6-1 in db versus phase coordinates. Figure 6-1 also 
contains the -1/N curves of Figure 5-1 for y = 1,38 x 10^, 

1.38 X 10® and 1.38 x 10^ 

With Y = 1.38 X 10’, the -1/N curve intersects the G(s) curves 
of the two systems at two points each. Among these the stable points 
for sustained oscillations are the ones on the left at the higher 
frequencies. The approximate magnitudes and frequencies of the 
oscillations are 6 x 10 ® rad and 4.4 rad/sec, respectively, for 
system 1, and 2 x 10 ® rad and 5.6 rad/sec, respectively, for system 
2. The curves in Figure 6.1 also show that for y considerably smaller 
than 1.38 x 10’, both systems will exhibit a stable response, although 
for certain values of y system 2 will show sustained oscillations 
while system 1 is stable. 

For the computer simulation, the input to the 1ST system, 
is set to zero, along with all the initial states, except for the 
vehicle position The initial value of 0^ is set at 5 x 10"® 
rad, which is chosen so that the input signal to the nonlinearity, 

6g, would be large enough to cause the torque to saturate, while 
at the same time the limiting value of the input signal is not 
exceeded. 

The following quantities are plotted from the simulation runs: 

Qy = vehicle position (rad) 
ijdy = vehicle velocity (rad/sec) 

= Gimbal position (rad) 
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Figure 6-1 








(jOq = Gimbal velocity (rad/sec) 

Tgp = Torque output of the nonlinearity (ft-lb) 

Error = Error input command (rad/sec) to the CMG 
= X - 

Figures 6-2 and 6-3 show the plots of the above listed quantities 
for system 1 with y = 1.38 x 10^. It may be noted from the plot of 
in Figure 6-3 that the system has a sustained oscillation. This 
oscillation is not seen on the other plots because of the large initial 
transients. Figure 6-4 through 6-5 show the continuation of Figures 6-2 
and 6-3 with proper scales. Figures 6-6 and 6-7 show the response plots 
for system 2 with y = 1.38 x 10^, and Figures 6-8 and 6-9 show the 
continuation of these plots with proper scales. The frequencies and 
magnitudes of oscillations obtained with the two systems are quite 
close to the predicted values. The small discrepancy is attributed to 
the discretization of the nonlinearity implementation on the digital 
computer. 

Figures 6-10 and 6-11 show the response plots for system 1 with 
y= 0.69 X 10^ and their continuations are shown in Figures 6-12 and 
6-13, respectively. Figures 6-14 and 6-15 show the response plots 
for system 1 with y = 1.38 x 10^ . As predicted, the system is stable 
for both of the lower y values. 
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7. Transfer Functions of the Sampled-Data LST System 


Since the actual LST system has sample-and-hold between the vehicle 
controller and the CMG controller, the system should be modelled as a 
sampled-data control system. Figure 7-1 shows the block diagram of the 
simplified LST system with sampled data. Since it is necessary to 
isolate the CMG nonlinearity from the linear dynamics for analytical 
purposes, a sample-and-hold Is inserted in front of the nonlinearity 
as an approximation. 

Referring to Figure 7-1, the following equations are written 
using e*, 0^* and Qy* as outputs. 


So‘^l‘^6'^7 


fi * - 




_Jl9 §. Q * ^ 

sAq J^g 







e* - N* 


. '^0 



(7-3) 


where the symbol * denotes the z-transform operation, and H* represents 
the discrete describing function of the CMG frictional nonli neari ty ; 

Equations (7-1) through (7-3) are portrayed by the sampled signal flow 
graph of Figure 7-2. Applying Mason's gain formula to this flow graph 
yields the determinant of the graph as 




Figure 7-1. A block diagram of the simplified LSI control system with sampled data 


OJ 











Figure 7“2. The Sampled Signal Flow Graph for the Equivalent System. 
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A = 1 + N* 


+ N* 
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The last equation is put into the form of (with N(z) = N*) 
1 + N(z)G(2) = 0 


where 


G(z ) - 


^ A^(z)-A2(z)A3(z)+A^(2)A-j(z) 
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-1 
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-1 
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Substitution of the system transfer functions into the above expressions 


yields 


JpS +K s+K, 

Aq = 
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-1 


^(z) = 0-z 
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,0pS^+K s^+K,Sj 

{ h P I J 


A2(z) = (l-z~^).,[ H(Kq+K^s) ) 


J^s^(JgS^+KpS+Kj) 
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A4(z) = (1-z ) J 


[ KjH(Kq+K^s) 
Jys3(jgs2+KpS+Kj)j 


(7-15) 


(7-16) 


The following system parameters are used for System, 1: 

H = 600 ft- lb-sec 
= 2.1 ft-lb-sec^ 

Kp = 5.75835 x 10^ 

K.| = 1.37102 X 10^ 

Kp = 216 ft-lb/rad/sec 
Kj = 9700 ft-lb/rad 
= 10^ ft-lb-sec^ 

Taking the z-transforms of the functions inside the brackets in 
Eqs. (7-13) through (7-16), we have the following results; 


A^(z) = + A 


14 


aT ^ '^15 
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2 -e 
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(7-17) 


a = - 51.429993 - j44. 42923 
a = - 51.429993 + j44. 42923 
A^ 3 = 1.0309376 x 10"'' 

= - 5.15469 X 10'^ - j5. 9669183 x 10'® 
= -5.15469 X lO'® + jS. 9669183 x 10'® 


Ayz) = A22 -y + A23 + 


24 


aT ^ ^'25 . ^ 


2 -e 


2 -e 


(7-18) 
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A ^2 = 3.5616776 x 10 
A 23 » 7.6870434 x 10''* 

A 2 ^ = -3.8435217 x 10'^ - j4. 8499764 x 10'^ 
Ajj - 3.8435217 x lO'^ + j4. 8499764 x 10'^ 


^3^^^ ^32 z-l‘ 


^33 


'34 


z-e 


Il-Jl -}- fl 
aT ^35 


z-1 


Z-G 


tr 


(7-19) 


A32 - 1.0 


^33" 


'34 


'35 


-2 


-2.2268891 x 10 
1.1134446 X 10"^ + jl. 6350537 x 10'^ 
1.1134446 X 10"^ - jl. 6350537 x 10'^ 


A4(2) 


Ml 


^M2 " 
^43" 


= A 


2(z-l)^ 


34.54837 

7.4564409 

■0.17352653 


- T .. ^ ^ JL- + + f\ 4- A 

^42 z-1 ^ ^43 ^44 ^ aT ^45 


z-e 


z-e^ 

(7-20) 


-2 


A.. = 8.6763263 x 10 ^ + jl. 6520832 x 10 


.-2 


M4 

^45 


8.6763263 x 10 ^ - jl. 6520832 x 10'^ 


The following system parameters are used for System 2. The 
same expressions for A^(z), A 2 (z), A 3 (z) and A 4 (z) are preserved. 
H = 200 ft-lb-sec 
Og = 3.7 ft-lb-sec^ 

Kq = 2 X 10^ 

= 3 X 10^ 
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Kp = 280 ft-lb/rad/sec 

Kj = 10 ^ ft-lb/rad 

5 2 

= 10 ^ ft-lb-sec 

The corresponding coefficients in Eqs. (7-17) through (7-20) 
are 

a = -37.83783 - j35. 651077 
a = -37.83783 + j35. 651077 
A ^3 = 9.9999976 x lO"^ 

-4:9999973 x 10"^ - j5. 3066848 x 10"^ 

A^ 5 = -4.9999973 x 10'® + j5. 3066848 x 10"^ 

A 22 = 3.9999932 x 10'^ 

A 23 = 4.8799929 x 10''^ 

Ap. = -2.43997 x 10'^ - j3. 150655 x 10"^ 

A 25 = -2.43997 X 10“^ + j3. 150655 x lO"'^ 

A 32 = 9.9999982 x lO"^ 

A 33 = -2.8 X 10"^ 

A 3 . = 1.4000002 X 10"^ + j8. 3391555 x 10"^ 

A 33 = 1.4000002 X 10'^ - j8. 3391555 x 10"^ 

A., = 39.999985 
A . 2 = 4.8800011 
A ^3 = -1.5144014 X 10'^ 

A^^ = 7.5720072 x lO'^ + jl. 1923421 x 10'^ 

A ^5 = 7.5720072 x lO"^ - jl. 1923421 x lO"^ 

It can be shown that if T approaches zero, the z-transfer function 
G(z) in Eq. (7-7) reverts to that of the continuous transfer function 
G(s) of Eq. (1-16). 



8. The Discrete Describing Function of the CMG Frictional Nonlinearity 

s 

In order to study the condition of self-sustained oscillations 
of the LST system with sampled data, it is necessary to evaluate 
the discrete describing function of the CMG frictional nonlinearity , 
N(z). 

The first step in the derivation of N(z) involves the inter- 
changing of the positions of the nonlinearity and the zero-order 
hold in Figure 8-la. This step is justified since the nonlinearity 
is amplitude dependent only, so that the signal of T^p is not affected 
by this interchange. Figure 8 -lb illustrates the transposition 
between NL and zoh. 

The second step involves the assumntion that 0^ is sinusoidal; 
that is , 

0g(t) = Acos(wt + (p) . (8-1) 


where A, w, and (p denote the amplitude, the frequency in radians, 
and the phase in degrees of the sinusoid, respecti vely. 

The z-transforra of Qg(t) is 


0g(z) = J Acos(-~- + (j))z 
k” 0 


-k 


(8-2) 


or in closed form 


Oq(z) = 


Az[(z - cos~-)cos(j) - sin^^intj)] 
z^ - 2zcos-^- + 1 


(8-3) 




(a) 



(b) 

Figure 8-1 

An important consideration is that because of the periodic 
nature of the sampler, 6g(t), 0g*(t), and Tgp*{t) are all periodic 
functions of period NT, where N is a positive integer ^2. 
Therefore, w = 2 tt/NT, and wT = 2 tt/N. 






The output of the nonlinearity in Figure 8 -lb is denoted by 
Tgp*(t) , and its z-transform is The discrete describing 

function (DDF) of the nonlinearity is defined as 


N(z) 


■ 83(21 


(8-4) 


It turns out that the discrete describing function (DDF) for 
N ~ 2 must be derived separately^ and a general expression for N(z) 
can be obtained for all N > 3. 


The DDF for N = 2 

Let Tgp(kT) denote the value of T^p*(t) at t = kT. For N = 2, 
the signal is a periodic function with a period of 2T. The 

z-transform of Tgp*(t) is written 

Tgp(z) = Tqp(0)(1 + + ...) + V^T)(z'^ + z‘^ + ...) 


TgF(0)z + Tgp{T)z 


(8-5) 


For the CMG frictional nonlinearity, it has been established 
in chapter 2 that 


^GF^^' ’’’gFO' 


p^ . 4 . 1 ~ |{1 - cos (tot + (j))] 
R~+ f * cos(u)t + 1 ) 1 )] 


0g<O (8-6) 



TgfrCt ) = TgpQ- 


|f^ " COS(wt + (J))] 
_ J .-^.. + _ COs(wt + jf ))] 


Let lis introduce the follov/ing notation: 

Tgp ( kT ) = ®G - ° 

V(kT) =Tgp(t)|t=kT ®g2-° 


We have 


V(kT) 


= T 


R + 1 


- |{1 “ COS(^^ + (j))] 


GFO 1 


R + 1 2^‘ 


cos{^+ 0)] 


'•'gf "^GFO' 


•r"TT + - cos(^+ (|))] 

+ |{] _ cos(^~+ ({))] 


For N = 2, 


TgpCO) Tgp (0) 

= T,p^(0) 

Tgp(T) = V(T) 
= V ( T ) 


0 <_ (|) < TT 


"IT < (f) < 2 tt 


TT ;< (J) < 27T 


0 ;< c}> < TT 


> 0 (8-7) 

(8-8) 

(8-9) 

<0 ( 8 - 10 ) 
>0 (8-11) 

( 8 - 12 ) 



Substituting U - 2 into Eq. (8-3), we have 


6e(z) = ^ 


(8-13) 


Using Eqs. (8-5) and (8-13), the DDF for N = 2 is determined. 


Tgp(0)2 + T„^(T) 




(8-14) 


Also, for N = 2, z = -1 , the last equation becomes 


N(z) = 


_ Tgp(O) - Tgp(T) 


2AcoS(|) 


(8-15) 


For stability analysis, we define 


F{z) = - 


1 _ 2Acos^ 

wTiT" Tgp(TT - Tgp(^ 


(8-16) 


The DDF for N > 3 

In general, the z-transform of the output of the nonlinearity may 
be written as 

00 M-i 

1 I Tgp(kT)z-'-"'N 

m=0 k=0 


N-1 


N-k 



(8-17) 



Using Eq, (8-2) for BgCz), the discrete describing function 
N(z) is written 


N(z) = 


T,3p(2) 


K, , 

I T3,(kT)z«-^ 


k=0 


(z^ “ "f) ? Acos(^ + (^)z 
k=0 ^ 


-k 


(8-18) 


The denominator of N(z) may be simplified as follows: 




k=0 

M-1 


= aIcos(?^^^)zN'^'1 


k=0 


(8-19) 


Thus 


I Tgp(kT)z'^'''-l 


N(z) = 


k=0 


aTcos(?^.1,)z«-^-1 

k=0 


(N > 3) 


( 8 - 20 ) 


As an alternative we may expand z^ - 1 as 


N-1 


z*' - 1 = H (z - 
k=0 


( 8 - 21 ) 


Then, using Eq. (8-3) for ©^(z), N(z) is written 



(8-22) 


N(z) = 


YT(5p(kT)z^'^[z2 _ 2zcos^+ 1] 

k<Ll 1 

n (z - e'^*^’^^^^)Az[(z - cos|^)cos(}) - sin|^ 1 n 4 )] 
k-0 


For N = 3, N(z) is simplified to 

. T3p(T)z . T3^(2T) 

J A(z ^l)X(z + 0.5)cos(|) - 0.866siri(|)] 


(8-23) 


For N > 3, 

^iVgp(kT)z^'^''' 

N(z) = ]^T 2 (^-24) 

A(z - 1) n (z - [(z - cos|^)coS({) - sin^in^)] 

k-2 

where in general, 

Tgp(kT) = V(kT) 

= Tgp^kT) 

where 0 < (27rk/N + c|)) £ 2 tt must be satisfied by appropriate conversion 
of the angle 2Tik/N + 

For stability studies the critical regions of F(z) = -1/N(z) 
should be constructed for N = 2, 3, with (^) varied from 0° to 
360°, and A from 0 to infinity. 


n ^ ^Ttk ... 

0 < + (f) < TT 

Tt < + (|) < 271 


(8-25) 



The following theorems on the properties of -1/N(z) are useful 
for simplifying the task of the construction of the critical regions. 

Tke.oA.ejn S-1 

To A any lnte.QAaZ H, the. maQyuXu.de. and pka^e. oi -J/W{z) Ae.pe.at 
ioA zoeAjy (j> = 2n/f>J AacUan6. 

Proof: The negative inverse of the discrete describing function 

is written 


F(z) = - 


-A[(z “ cosl^)coS(l) - sin|^int}3](z^ - 1) 


NTiT 


N-1 

I Tgp(kT)z 


N-k~l 


k=0 


(z^ - 2zcos^ + 1 ) 


(8-26) 


Let 


FJz)Fp(z) 
" F3(2)F^(2) 


(8-27) 


where 


F.j(z) = -A[(z - cos^)cos(f) - sin^in<|)] (8-28) 

T^{z] = z^'^ - 1 (8-29) 

T.iz) = I Tpp(kT)z^'*'”^ (8-30) 

k=0 

F^(z) ~ {z - 2zcos^— + ]) (8-31) 
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Let 


F(z)„ = [F(z)I^ ^ 


Fi(z)n = [F^Cz)]^ ^ ^ ^ 2,r/N 
= (p + 2 tt/N 


Then, 



Fi(z) 

u 

F3(z)j^ 


p 



and 


Arg[F(z)] - Arg[F(z)j^] = Arg[F^(z)] - Arg[F3(z)] 

- Arg[F^(z)j^] + Arg[F3(z)j^ 


Also, 


F^(z) 


Asinl^ 




Asin~: = |F^(z) 


Let us express F,(z)jj as 


(8-32) 
(8-33) 
(8-34) 

(8-35) 

] (8-36) 

(8-37) 
(8-38) 
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(8-39) 

where 

= [Tgp(kT)]^ ^ ^ 2TT/N (8-40) 

It can be shown that for any integral N, 

N-1 

F3(z)^ = zF3(z) = I Tgp(kT)z'^-^ (8-41) 

Ic ™ 0 


Then 5 Eq. (8-35) becomes 



zFgCz) 


= 1 


(8-42) 


The argument of F^(z) is 

Arg[F-|(2)] = ({) - tt/2 


(8-43) 


Then 


Arg[F-j(z)j^] = (|) - tt/2 + 2tt/N 


(8-44) 


Thus, Eq. (8-36) becomes 



Ar^[F{z)] - Arg[F(z)^] ^ - Arg[F 3 (z)] + Arg[F 2 (z)pj] 

= - ~ - Arg[F3(z)] + ArgCzF^Cz)] 

= - “ Arg[F3(z)] + Arg[F3(z)] 

= 0 

Q.E.D. 

As an illustrative example of Theorem 8-1, let us consider the 
case of N = 4. 

Let I = yTgpoA = K 

Then 


T -n,Tl - R + K (R + l)[cosU_+_2irk/N) - Mj 
'GF ■ 1 - kTr +~r)Tcos($ + 2 irk/N) - TT GFO 


Tgp"(kT) = 


R - K(R - 1 

)rcos(.j) + 2irk/N) 

- 

1 - K(R - 1 

)[cos (J) + 2Trk/N] 



For 0 < <j) < IT, 


F3(z) = Tgp'(0)z^ + + T^p+OT) 

" ^GF ^GF "''gF ^GF 


(8-45) 

(8-46) 

(8-47) 

(8-48) 

(8-49) 


It is easy to see that 
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^GF ''’gF 

Tgf'(T)n = T3p"(2T) 
Tgf"( 2T), = Tgp"(3T) 
%"(3T), - T3p-(0) 


which proves that ^ 2 ( 2 ) 1 ^ zF 2 (z). 

Similar results are obtained for tt < <|) < 2ir. 

Tke.o^<im S-2 

Fo/l odd iV fW ^ 3 ), ddie. niagyuXude and pfL06e, 0^ - 7 /A/fzJ ^e.pe.a-t 

^ 0 ^ c.veAj/ c|> = 7 t/M. 

Proof: Let f^(z)jyj> defined as F(z) , 

F-j(z), and F^tz) with cp replaced by (f> + tt/N, respectively. 
Then , 


Arg[FT(2)fj] = <(> - | + | 

Ar-g[F(z)] - Arg[F(z)j^] = - ^ + Arg[F 3 (z)j^] - Arg[F 3 (z)] 

Using the same notation as in Theorem 8-1, it can be shown that 
for odd N > 3, 


(8-50) 


(8-51) 

(8-52) 



(8-53) 


F3(2)n = 


-(N-D/2 


F3(2) 


Thus, 


Arg[F3(z)^] = iT - Arg[F3(z)] 

= I + Arg[F3(z)] 


(8-54) 


Again, 



(8-55) 


and 


Arg[F(z)] - Arg[F(z)jj] = - + Arg[F 3 (z)j^] - Arg[F 3 (z)] 

= 0 (8-56) 

Q.E.D. 

As an illustrative example of Theorem 8-2, consider the case N = 3. 

For 0 < (|> < 7T, 

FsCz) = Tgp‘(0)z2 + Tgp-(T)z + Tgp+(2T) (8-57) 

+ Tsp"(2T)^ (8-58) 


It can readily be shown that 
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"'^GF "" ■''^GF 

"''gf " ''''gf 

'^GF " '^GF 


(8-59) 


To carry out one of the identities above. 


_ _ R - K( R “ 1 ^fcos (0 + tt) - 

*GF '^‘^N GFO 1 - “K{R lT[cosTr+ 

T ~{r\\ = r R + K ( R + 1)(cos<t - 1) 

'gf GFO 1 - K(nn)(cos^ - 1) 


R + K(R - 1)(cosj) + 1). 
1 + K(R - l)(cosii + 1) ' 


GFO 
(8-60) 

(8-61) 


Using the relation, R + KR - 1 = 0, we have 


^GF “ ^GF 


(8-62) 


Thus , 




(8-63) 


The significance of the last two theorems is that the critical 
regions of -'1/N(z) need be computed only for 0^ < c|) ;< tt/N for odd 
N, and 0° < <_ 2 tt/N for even N. 



51 


S'- 5. 

Aiympotlc, BtkavloK ~l/hJ[z) <x6 A appA.oadieJ> l¥i{^nLty, 

(a) ^n?|-l/N(z) I = » 

A-x» 

(b) For even N ^ 4, 0 ^ 2 tt/N, 

Urn Arg[-1/N(z)] = (~ - i)ir + (}> 

A->oo 

For odd N > 3, 0 < (j> < tt/N. 
tun Arg[-1/N(z)] = (1 - + 4> 

A-Hxj 

(c) For N=2 0£(f)£Tr 

tim Arg[-1/N(z)] = 0° 0 5 < ‘n-/2 

A-^ 

= 7T TT/2 ^ < 7T 

(d) For N = 3 0^4^^ 'it/3 

Urn Arg[-1/N(z)] = - ~ + (|) 

A-h» 

Proof: 

We can easily show that 

Urn Tgp (kT) = -TgpQ ^ 

A-«o 

lim Tgp'^ 0<T) = TgpQ 6 >0 


(8-64) 


(8-65) 


(8-66) 


(8-67) 


( 8 - 68 ) 


(8-69) 


(8-70) 



The magnitude of -VN(z) is directly proportional to A as A 


approaches infinity, thus (a) is proved. 
For N = 2 0 < 4> < tt/2 


■ -'IXM = 


GF 


Thus, 


tim Arg[F(z)] = Urn Arg L - — HTT 

A-xx> A-^ 


= Arg =0° 0 < (}) < 7 t/2 

'gfo 

~ ■n' 7r/2 £ 4> < 


This proves, item (c). 
For N = 3 0 < (|) < 7 t/3 


P( 2 ) = -A[(z + 0.5)cos({) - 0,866sin(j)1(z - 

T^p(0)z2 + Tgp(T)z + T^_^{2T) 


GF' 


tim Arg[F(z)] - 

/\-w 


5it 9 

T " ^ Arg[Tgp(0)z + Tgp(T)z + Tgp(2T)] + 


1^ - Arg[-z^ - z + 1] + 4 > 


A A 

y-+ 


(8-71) 


(8-72) 


(8-73) 


(8-74) 


This proves item (d). 
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For N ^ 4 and even 0 < < Ztt/N 

Using Eqs. (8-27) through (8-31), we have 

tern Arg[F(z)] = Urn [Arg[F,(z)] + Arg[Fp(z)J - Arg[F.(z)] - Arg[F.(z)] 
A-^ A-^ ^ ^ 



tlm ^Arg[F3(z)] 


+ Arg[Fp(z)] - Arg[F.(z)] 

(8-75) 


Arg[F 2 (z)] - Arg[F,(z)] = Arg[(z - 1) n (z - 

^ k=2 

= Arg[e'^^'^^(e'^*'^'^^ - e"'^*^^^) H (8-76) 

k=2 

Arg[p 2 (z)] - Arg[F,(z)] = H (1 - 

k=l 

. (N . 3)(|I^) + 

^ ^=1 

- - 3)(|^) + (N - 4)(- f) + (8_77) 


llm Arg[F^(z)] = tor] Arg ^ Tpr(kT)z^’’^’’^ 
A-^ k=0 
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N/2-1 M 1. 1 N-1 , 

Arg I - TgpoZ^“ “ + Arg J 


k=0 


k=N/2 


GFO' 


N/2-1 


Arg I - 


(t- 1)^ 


k=0 


GFO 


TT TT 

2 ' N 


Thus, 


^ Arg[F(z)] = $ - I + 1+ (N - 4)(- J) + (N - 3)(^) + 


- n N A ^ k, ^ ^ 
2 ' N 


(2 - |l)’f + <l> 


For N > 3 and odd 0 5 4> 5 'tt/N 


For this case, 


Lim 

A-^ 


(N-D/2 


im Arg[F3(z)] = Arg J + Arg’^’ T 


N-1 


k-0 


GFO^ 


N-k-1 


k^i+))/Z 


_ (N - 3)7t 
2N 


(8-78) 


TT TT 

2^H 


(8-79) 


(8-80) 


Thus, 



^Arg[F(z)] = ^ - | + 1 + (N - 4)(- ^ + (N - 3)(|^) + V| 

w 7 N-3 N 

= (3 - g- - 2)f + j^)ir + it> 

= 0 - ^)f + <i> 


and (b) is proved. Q.E.D. 

Tk<Loh<m S-4. 

Aiyrr\p^otlc. 'Stkcvlofi o{^ -1 /U[z] A app^oachsJ) Z2Ao, 

■Urn F(z) Ly 

■ "^GFO 


for al 1 if and all N. 
Proof: From Eq. (8-20) , 


F(z) = - 


1 

NlTT 


-AVcos(-2fL.,),N-k-l 


k=0 


N-1 . 

I Tgp(kT)z^ 

k=0 


l-k-1 


Lim F(z) = 
/HO 


'PU 

- l cos(~^ + (f))z 
k=0 

N-1 T„(kT) 
Urn I 
A-+0 k=0 


A 


N-k-1 


M-k-1 


(N - 3)tt 
2N 


(8-81) 


(8-82) 


(8-83) 



(8-84) 


N-1 

- I cos{?^+ <^)z 


k=0 


1-1 T.7rkT) 

I Um z 


k=0 A-^0 


A 


N-k-1 


N-k-l 


T. (kT) 

Therefore the -problem is - th-at of finding ^ 

A->0 ^ 

First, let Tgp(kT) = T^p"(kT). Then 


tim 

A^O 


Tgp'(kT) 

A 


T A(R + IT ~ ^^GFO 

A-0 YTgpoA - YTgpoAcos(^ + 


" '''gFO A(R + 1) ^ 


GFO 


[cos( ^ 


+ 4^) - 1] 


where the fact that tlm [1/(R + 1)] = 1 has been used. 

A^O 




^ = Zbn Y 
A->0 ^ A-^ ^ 


/Ay^T^FO^^^ - 1 


2yTqfo^ 


4y^T 


GFO 


2 2 
A 


Or, 


aTFTT) 


Ujm -T 
A-^0 ^ 


1 


2yTgFo^ 


^’'''’'gfo'^ 


1 j. /I 2-r 2.2 

1 + 4Y Tgpg A 


(8-85) 


(8-86) 


J 


(8-87) 
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Expandi ng 



2 2 2 


into a power series, and using only the 


-f-i rst-two terms., we- have., 


JU.m 

A-vO 


+ Ij 


/I 2^ 2„2 

1 i Tppn ^ 

= Um i 


A-K) 


A| 




'^'^GFO 


(8-88) 


Thus , 


Zlm 

A^O 


Tgp-(kT) 

A 


+ ^YF0^tcos(^+ ^) - 1] 


_ y 2 /277k , 

yTqfo N 


(8-89) 


Similarly, it can be shown that 


tim 

A"^0 


YphkT) 

A 


- Ztm 

A-^0 


TGp“(kT) 

A 


^ T 2 /2TTk . . \ 

^Yfo 


Now , 


tdn F(z) = 
A"^0 


N-1 

Y / 277k 
" I cos(^ + (j)}z 

k=0 


,N-1 


■ 2Trk 




k=0 


1 



(8-90) 


(8-91) 


for all (j) and all N. Q.E.D. 



